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^ ■ Abstract 

^ I We consider the processes of compare the heavy quark production 

using the unintegrated gluon distributions. The numerical predictions 
'nI" , for high energy nucleon-nucleon and photon-nucleon cohisions of the 

fcr-factorization approach (semihard theory) are compared with the ex- 
^SJ ' perimental data from Tevatron-colhder and HERA. The total produc- 

^" . tion cross sections and px distributions are considered and they are in 

ly-v \ reasonable agreement with the data for rather value of QCD scale. 
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1 Introduction 

The investigation of heavy quark production in high energy hadron colhsions is 
an important method for studying the quark-gluon structure of hadrons. The 
description of hard interactions in hadron colhsions within the framework of 
QCD is possible only with the help of some phenomenology, which reduces the 
hadron-hadron interaction to the parton-parton one via the formalism of the 
hadron structure functions. The cross sections of hard processes in hadron- 
hadron interactions can be written as the convolutions squared matrix elements 
of the sub-process calculated within the framework of QCD, with the parton 
distributions in the colliding hadrons. 

The most popular and technically simplest approach is the so-called QCD 
coUinear approximation, or parton model (PM). In this model all particles in- 
volved are assumed to be on mass shell, carrying only longitudinal momenta, 
and the cross section is averaged over two transverse polarizations of the in- 
cident gluons. The virtualities q^ of the initial partons are taken into account 
only through their structure functions. The cross sections of QCD subpro- 
cess are calculated usually in the next to leading order (NLO) of as series 
[1, 2, 3, 4]. The transverse momenta of the incident partons are neglected 
in the QCD matrix elements. This is the direct analogy of the Weizsaecker- 
Williams approximation in QED. 

Another possibility to incorporate the incident parton transverse momenta 
is referred to as /c^-factorization approach [5, 6, 7, 8, 9, 10], or the theory 
of semihard interactions [11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21]. Here the 
Feynman diagrams are calculated taking account of the virtualities and of all 
possible polarizations of the incident partons. In the small x domain there 
are no grounds to neglect the transverse momenta of the gluons, qix and g2T, 
in comparison with the quark mass and transverse momenta, pix- Moreover, 
at very high energies and very high piT the main contribution to the cross 
sections comes from the region of qiT ~ Pit or g2T ~ Pit, see [22, 23, 24] for 
details. The QCD matrix elements of the sub-processes are rather complicated 
in such an approach. We have calculated them in the LO. On the other hand, 
the multiple emission of soft gluons is included here. That is why the question 
arises as to which approach is more constructive. 

In Sect. 2 we present the cross sections of heavy quark production in 
hadron-nucleon and photon-nucleon collisions in the /c^-factorization approach. 
This approach accounts for the diagrams without neglecting by the virtuali- 
ties of incident partons (mainly gluons), their polarization, etc. The detailed 



discussion of the /cy-factorization approach, and its numerical difference from 
the parton model can be found in [23, 24]. 

The different approaches for the unintegrated gluon distributions are dis- 
cussed in Sect. 3 and 4. We show that the difference between them is not very 
large, and after integration the unintegrated gluon distributions reproduce the 
initial structure functions. At the moment we have no more realistic unin- 
tegrated parton distributions which are extracted from the data by the same 
way as structure functions, i.e. by using the evolution equations, boundary 
conditions at some small Qq, etc. 

The experimental data on heavy quark production at Tevatron-collider 
and HERA are discussed in Sect. 5. Our predictions with the selected kx- 
factorization (which seens to be rather natural) are in reasonable agreement 
with pT-distributions of 6-quarks produced in FNAL, as well as with the data 
on total cross sections of charm and beauty production at HERA. Predictions 
of the /cT-factorization approach for heavy quark production at LHC are also 
given. 

2 Heavy quark production in hadron-hadron 
and photon-hadron collisions in /cy- factorization 
approach 

The conventional parton model expression for the calculation of heavy quark 
hadroproduction cross sections has the factorized form [25]: 

a{ab -^ QQ) = Y1 j dxidxjga/i{xi, HF)9b/j{xj, fiF)o-{ij -^ QQ) , (1) 

where ga/i{xi, fip) and gb/j{xj, fip) are the structure functions of partons i and j 
in the colliding hadrons a and b, hf is the factorization scale (i.e. virtualities of 
incident partons) and da{ij -^ QQ) is the cross section of the direct subprocess 
{ij -^ QQ which is calculated in perturbative QCD. The last cross section can 
be written as a sum of LO and NLO contributions, and its analitical form can 
be found in [1, 2, 3, 4]. 

The principal problem of parton model is the collinear approximation. The 
transverse momenta of the incident partons, Qix and Qjt are assumed to be zero, 
and their virtualities are accounted through the structure functions only; the 



cross section a{ij -^ QQ) is assumed to be independent on these virtualities. 
Naturally, this approximation essentially simplifies calculations. 

Sometimes namely these simplifications result to the serious disagreement 
with the experimental data. 

Without the discussed simplifications of PM, the differential cross section 
of heavy quark hadroproduction in LO QCD (Fig. 1) has the following form 
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Here we have used the Sudakov (light cone) decomposition for quark mo- 
menta pi^2 through the momenta of colliding hadrons pa and Pb-, 
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(The more accurate relations are q\ = — g^j./(l — y), g| = — Q'2t/(1 ^ ^) but 
we are working in the kinematics where x,y << 1). 

qi 2T are the gluon transverse momenta; 



xi = —^e ^i , yi = ^F=^e^i, 



m2T _„. rn2T „* ,„s 

X2 = -y=^e y^ , 1/2 = -/^e^'' (6) 



2 _ 2,2 

^^1,2T — Tn -\- Pi2Ti 



^We put the argument of as to be equal to gluon virtual! ty, which is very close to the 
BLM scheme[26]; (see also [14]). 



where yl2 ^-^e the quark rapidities in the hadron-hadron c.m.s. frame and m 
is the mass of produced heavy quark and we take these masses 

me = 1.4 GeV , rrn, = 4.6 GeV , (7) 

for the values of short-distance perturbative quark masses [27, 28]. |MqqP is 
the square of the matrix element for the heavy quark pair hadroproduction. It 
is calculated in the Born approximation of QCD without standard simplifica- 
tions of the parton model. Contrary to the mention of [18], the transformation 
Jacobian from x,y to yl,y2 is accounted in our matrix element. 

The unintegrated gluon distributions f g{y , qir , jj) and fg{x,q2T, fJ') will be 
discussed in details in the Sect. 3. 

To provide the kr factorization of eq. (2) the cross-section of the elementary 
subprocess for the heavy quarks production, gg — > QQ was calculated in the 
axial gauge where the spin part of the gluon (gj) propagator takes the form 



d^,u{qi) = Qfiu 






with the gauge vector n^, = (pi +P2)u — (Pi +P2)t,u (that is rii, = (no, 0, 0, 0) in 
the frame where the longitudinal part of quark pair momentum {pi +^2)2 = 0). 

In this gauge both the DGLAP and BFKL leading logarithms in the un- 
integrated distributions fg are given by the ladder type diagrams shown in 
Fig.l without any extra interactions between the upper and lower parts of the 
graphs (that is between the functions fg{x^q2T, jj) and fg{y,qiT, f^))- 

Being averaged over the colours of the incomimg gluons the matrix element 
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Here iVc = 3 is the number of colors, the variables s, t, -u are defined for 
99 ~^ QQ subprocess. 
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The elements rriik are the numerators of the lowest order QCD Born diagrams 
squared and averaged over quark helicities; the incoming transverse momenta 
gi 2T play the role of gluon polarization vectors, 
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V3 = qiT ■ q2T (gi - 92) - ((2gi + ^2) ■ q2T) qit + ((2^2 + gi) ■ qir) q2T 

is the contribution of a triple gluon vertex. 

Thus we can calculate the heavy quark production cross section without 
standart simplifications of the parton model since in the small x region there 
are no grounds for neglecting the transverse momenta of the gluons gi,2T in 
comparison with the quark mass and transverse momenta J9i,2t- Although 
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explicit expressions for rriik are rather bulky, they can be easily computed us- 
ing any analytical computer program (say, REDUCE or MATHEMATICA). 
For example, the dependence of the matrix elements on the gluon transverse 
momenta is a main source for the nontrivial azimuthal correlations arising be- 
tween quarks momenta pi,2T which could not be described in the conventional 
parton model. 

In the case of high energy heavy quark photoproduction on the proton 
target we should consider two possibilities, the same resolved production via 
quark-gluon structure of the photon. Fig. 1, where one proton should by re- 
placed by photon (with its parton distribution, different from the case of in- 
cident nucleon), and direct production in the photon-gluon fusion, 7p -^ QQ, 
Fig. 2, where the incident photon has the fixed momentum. 

The cross section of the direct interaction can be written by the same way 
as Eq. (2) where the unintegrated gluon distribution fg{y,qiT, fJ') should be 
replaced by S{y — 1). The cross section of the direct interaction can be written 
as 



d'^PiT 71 J qf 



2^ 



. ^l(i-.)^..^(f.^).™^(i;-^)^(0) 

Di=piT + mQ, D2 = {qr - Pirf + rriQ , (10) 

here aem = 1/137 and cq is the electric charge of the produced heavy quark. 



3 Unintegrated gluon distributions in differ- 
ent approaches 

The unintegrated parton distribution fa{x,qT,fi) determines the probability 
to find a parton a initiating the hard process with the transverse momentum 
qrp (and with factorization scale /i). 

At very low x, that is to leading log(l/x) accuracy, for the case of gr of 
the order of scale /i the unintegrating gluon distribution are approximately 
determined [11] via the derivative of the usual structure function: 

f,{x,qr,f^) ^Ux.q') = ^if^- (11) 



However at a large x the unintegrated density (11) becomes negative. 



To restore the function fa{x,qT, f^) on the basis of the conventional (inte- 
grated) parton density a{x, A^) we have to consider the DGLAP evolution^ 
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Here, a{x, qf) denotes xg{x, qf) or xq{x, qf), and Paa' are the splitting functions. 

The first term on the right-hand side of Eq. (12) describes the number 
of partons 6a emitted in the interval A^ < g|. < A^ + (5A^, while the second 
(virtual) term refiects the fact that the parton a disappears after the splitting. 

The second contribution may be resummed to give the survival probability 
Tg that gluon g with the transverse momentum g^ remains untouched in the 
evolution up to the factorization scale 



Tg{qT,^) =exp 
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Thus, the unintegrated gluon distribution fg{x,qT,fi) has the form 



fg{x,qT,fi) =Y1 
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where the cut-off A is used in Eqs. (12-14) [32, 33]. 

The expression (14) with the survival probability (13) provides the positiv- 
ity of the unintegrated probability fg{x, qx, jj) in the whole interval < a; < 1. 
It is necessary to note that the calculation of the integrals in Eqs. (13) and 
(14) should be fulfilled with rather high and equal accuracy because the factor 
Tg{qT,fi) compensate the singularity at A ^ 1 in Eq. (14). 

For the case of one loop QCD running coupling ^^(pt) = 47r/(61np^/A^) 
the factor Ta{qT,f^) can be written down explicitely. In particular, the gluon 
survival probability (which enters our formulae) reads: 
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^For the g -^ gg splitting we need to insert a factor z' in the last integral of Eq. (12) to 
account for the identity of the produced gluons. 
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where 

Ek{qT,f^) = (-] [Et{k\n{fi/A)) - Et{k\n{qT/A))] , (16) 

and the integral exponent Ei{z) = Jl^ y expt; np and Nc are the number of 
light quark flavoures and the number of colours, respectively, and b = 11 — ^np- 

In the leading log(l/a;) (i.e. BFKL) limit the virtual DGLAP contribution 
is neglected. The scale // is of the order of qx- So T^ = 1 and one comes back 
to Eq. (11) 

/r^-(x,,^,;.) = ^^^ ,X = gT. (17) 

In the double log limit Eq. (12) can be written in the form 



d 
fa^^ix,qT,IJ,) = -gj^ [a(x,A2)T„(A,/i) 



A=(JT 



which was firstly proposed by [32]. In this limit the derivative dTa/ din X"^ 
cancels the second term of the r.h.s. of Eq. (12) (see [33] for a more detailed 
discussion)^. 

Finally, the probability fa{x, qr, 1^) is related to the BFKL function ip{x, g|.) 

by 

^(a:,g^) = 4y2 7rVa(x, gT,/u)/g| . (19) 

Note that due to a virtual DGLAP contribution the derivative da{x, \^) /d}? 
can be negative for not small enough x values. This shortcoming of Eq. (17) 
is partly overcome in the case of Eq. (18). 

We have to emphasize that fg{x,qT,fi) is just the quantity which enters 
into the Feynman diagrams. The distributions fg{x,qT, fi) involve two hard 

^ There is a cancellation between the real and virtual soft gluon DL contributions in the 
DGLAP equation, written for the integrated partons (including all fcy < /i). The emission of 
a soft gluon with momentum fraction (1 -- z) — ^ does not affect the x-distribution of parent 
partons. Thus the virtual and real contributions originated from 1/(1 — z) singularity of 
the splitting function P{z) cancel each other. On the contrary, in the unintegrated case the 
emission of soft gluon (with q'rp > kr) alters the transverse momentum of parent (i-channel) 
parton. 
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scales^: qt and the scale /x of the probe. The scale /i plays a dual role. On the 
one hand it acts as the factorization scale, while on the other hand it controls 
the angular ordering of the partons emitted in the evolution [34, 35, 36]. The 
factorization scale /i separates the partons associated with the emission off both 
the beam and target protons (in pp collisions) and off the hard subprocess. 
For example, it separates emissions off the beam (with polar angle 9 < 90° 
in c.s.m.) from those off the target (with 6 > 90" in c.s.m.), and from the 
intermediate partons from the hard subprocess. This separation was proved 
in [34, 35, 36] and originates from the destructive interference of the different 
emission amplitudes (Feynman diagrams) in the angular boundary regions. 

If the longitudinal momentum fraction is fixed by the hard subprocess, then 
the limits of the angles can be expresseed in terms of the factorization scale 
Ij, which corresponds to the upper limit of the allowed value of the s-channel 
parton kx- 



4 Numerical values of unintegrated gluon dis- 
tributions 

As was shown above, there exist several ways for estimation the gluon unin- 
tegrated distributions. Now we will compare numerically several of them to 
see the theoretical uncertainties which should result in the uncertainties in the 
predictions for cross sections of heavy quark production. 

First of all, it is evident that unintegrated gluon distributions obtained 
with the help of Eqs. (12)-(14) should coincide after their integration over g^ 
with the used structure functions. 



,q')=J%f9{x,qT,fi). (20) 



dq 
xgix "■" ' 

■ q 



Some problem is that the values of fg(x,qt,n) are principally unknown at 
very small g|, ~ g^. So instead of Eq. (20) we can write 

f do 
xg{x, q^) = xg{x, ql) + / -^fg{x, qr, /i) . (21) 

^This property is hidden in the conventional parton distributions as qt is integrated up 
to the scale /x. 
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with g2 ^ 1 GeV^. 

The calculated values of heavy quark production cross sections depend [38] 
on the analytical form of the cut-off parameter A in Eqs. (12)-(14). However, 
if Eqs. (13) and (14) are calculated with the same numerical accuracy, this 
dependence becomes rather weak. 

Now we compare the results for sum rules, Eq. (21), for two cases, 

A = — ^^ (22) 

fi + QT 



^ (23) 



and 

A = - 

V/^^ + Qt 

The results of calculations are presented in Fig. 3a and one can see that 
the variants (22) and (23) give practically the same results. In the futher 
numerical calculations we will use the variant (22) which corresponds to the 
angular ordering in the gluon emission [34, 35, 37]. 

In Fig. 3b we compare the q^ depencences of unintegrated gluon distribu- 
tions given by Eq.(14) [33], the KMS distribution [39], which was based on the 
BFKL equation and fitted to the F2 HERA data just in unintegrated form, 
and simplified formula (11) [11]. In the first case we present the results for the 
value of QCD scale /x^ in Eqs. (13), (14) equal to 100 GeV^ (i.e. about 4m^). 
Here we have used rather more or less realistic gluon distribution GRV94 HO 
[40] which is compatible with the most recent data, see discussion in Ref. [41]. 
For a low x all three distributions are rather close to each other while at a 
larger x ~ 0.05 the KMR prescription gives for g^ > 10 GeV^ about twice 
larger gluon density. 

Eq. (2) enables us to calculate straightforwardly all distributions concerning 
one-particle or pair production. One-particle calculations as well as correla- 
tions between two produced heavy quarks can be easily done using, say, the 
VEGAS code [42]. 

However there exists a principle problem coming from the infrared region. 
Since the functions ip{x,qf) as well as fg{x,q^^, fj,) are unknown at small val- 
ues of g| and qf, i.e. in nonperturbative domain we calculate separately the 
contributions from qf < Ql, ql < Q^, qf > Ql and qf > Ql [22, 23, 24]. 
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0.1 12 

where the unintegrated gluon distributions /c,(x, gy, jj) are taken from Eq. (14). 
In the numerical calculations we use the values /i^ = rn^, /i^ = Am\ and /i^ = s. 

The first contribution in Eq. (24) {ql < Q^, g| < Q^) with the matrix 
element averaged over directions of the two-dimensional vectors qiT and q2T is 
exactly the same as the conventional LO parton model expression. It is multi- 
plied by the 'survival' probability T'^{QQ,fi'^) not to have transverse momenta 
qiT, q2T > Qo- We assume Ql = 1 GeV^. The sum of the produced heavy 
quark momenta is taken to be exactly zero here. 

The next three terms (when one, or both q^^ > Ql) contain the correc- 
tions to the parton model due to account the different gluon polarizations, 
their virtualities and transverse momenta in the matrix element. The relative 
contribution of these corrections strongly depends on the initial energy. If it 
is not high enough, the first term dominates, and all results are similar to the 
conventional LO parton model predictions [23]. In the case of very high energy 
the opposite situation takes place, the first term can be considered as a small 
correction and our results differ from the conventional ones. So we need the 
highest energies for investigation the /c^-factorization effects. 
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5 Heavy quark production in the /ct- factorization 
approach 

Let us compare the results of our numerical calculations with the data on 
beauty production at Tevatron-collider. In Fig. 4a we present the data [43] 
on b quark production with pt > Pmin identified by its muon decay, as the 
function oi pmin at 1.8 TeV. The curves of different type show the calculated 
results with different scales /i^ in (13) and (14), /x^ = rri^, ij? = 4m^ and 
fi^ = s, respectively, where s is the invariant energy of the produced heavy 
quark pair, s = xys. The calculated values of one-particle Ex distributions, 
da/dEx, in /c-r-factorization approach are presented in Fig. 4b together with 
the data [44] extracted from muon tagged jets. In both Figs. 4a and 4b our 
results are slightly smaller than the data. 

The similar data on bb production with pt > Pmin obtained by UAl (circles) 
and CDF Coll. (Triangle point) at ^/s = 630 GeV taken from [45] are presented 
in Fig. 4c together with the results of our calculations. Again, our results 
are in agreement with UAl and slightly smaller than the CDF data. Some 
disagreement between UAl and CDF experimental data is evident. 

The experimental data for charm production [46] at 1.96 TeV are compared 
with our calculations in Fig. 5. Now our curves slightly overestimate the data. 



Table 1 The total cross sections of charm and beauty production in the 
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rrir 





all rapidities 
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cc 


bb 


cc 


bb 


14 TeV 


18.8 mb 


0.92 mb 


4.24 mb 


0.245 mb 


1.96 TeV 


2.95 mb 


89.6 /ib 


843 /ib 


33.3 /ib 


1.8 TeV 


2.7 mb 


84.5 ^ih 


780 /ib 


30.2 /ib 



The energy dependences of the total cross sections of cc and bb pair pro- 
duction are presented in Table 1. 

As it was mentioned, at comparatively small energies the first term in 
Eq. (24) dominates and the result of /cj-factorization approach practically 
coinsides with LO parton model prediction. With increasing of the energy the 
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difference between the considered approaches increases and at the LHC energy 
the /c^-factorization predictions are about two times larger than the LO parton 
model ones (i.e. of the order of NLO parton model predictions). Generally, 
the /cT-factorization yields a more strong energy dependences both for cc and 
hh production. It can be explained by additional contributions appearing at 
very high energies in the /cy-factorization approach, see [22]. The main part 
of these contributions corresponds to the configurations where the transverse 
momentum of heavy quark is balanced not by another heavy quark but by the 
qx of gluons. 

Note that we underestimate the 6-quark cross section but overestimate 
the charm production. This may be caused either by the behaviour of the 
gluon densities used in the calculations (too low density in the kinematical 
region (x ~ 0.01, /x^ ~ 50 GeV^) corresponding to the beauty production 
at the Tevatron and too high density at a lower scale and x relevant for the 
charm production) or this may indicate some experimental problem in selecting 
the 6-quark events. Not excluded that sometimes the muon or jet from a 
charm production (which has a much larger cross section) was treated as the 
lepton/jet coming from the beauty. 

It is not quite clear what the renormalization scale /i/j should be used 
as the argument of QCD coupling a^ in Eqs. (2) and (24) within the LO 
approximationof /c^-factorization approach. As the main version in the present 
paper we choose Hr = Q^, motivated by BLM prescription [26]. On the other 
hand, based on the explicit one- loop calculation in Ref. [32], the argument was 
given in favour of the running coupling depending on the largest virtuality in 
the ladder cell under consideration^. 

By this reason in Fig. 4a we present the curve corresponding to the /i/j = 
rriT (that is in Eq. (24) we put the coupling asimx))- For the case of beauty 
production this choice reduces the predicteed cross section of about 2-^3 times 
for pTmin = 5-30 GcV/c, respectively. However for a lighter charm quark using 
the asirriT) we obtain the cross section only 1.2 -^ 1.5 times smaller than that 
with as{qf). 

Heavy quark production at the Tevatron within the framework of kr- 
factorization was studied recently in [47]. The authors used the unintegrated 
gluon densities obtained from the Linked Dipole Chain model. With the " stan- 
dard" gluons they obtained the results close to our calculations and underesti- 



^Thc prcexponent factor lii(^/A)/ln(gT/A) in Eq. (15) effectively replaces the one-loop 
running coupling as{p) by the as^qx), close to the BLM prescription [26]. 
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mated the beauty cross section extracted from the measurement of the muon 
and jet distributions. To describe the data it was needed to take so-called 
"leading" gluons, that is the version of the model which neglects the light 
quark contribution and the non-singular terms in the splitting functions. 
We see no possibility to justify such an approach. 

On the other hand the new CDF data where the S-meson was identified 
via the B — > J/\l/ X decay mode give a lower cross section. The last analysis 
of this new CDF data gives 29.4l;5;|/ib [48] in agreement with our result (25/ib) 
for beauty production at 1.96 TeV for \yh\ < 1. 

The energy dependences of the total cross sections of cc and bb photopro- 
duction are presented in Fig. 6. The data on charm (Fig. 6a) and beauty 
(Fig. 6b) photoproduction cross sections are taken from [49, 50, 51]. All three 
values, used for the scale give very similar results which are again slightly 
lower than the high energy data. Here we show separately the sum of direct 
and resolved contributions, as well as the only resolved ones. The resolved con- 
tributions have more strong energy dependence. At energies W ~ 300 GeV 
they give 12 — 15% of the total cross sections but in some kinematical regions 
they can even dominate [38], for example, in the target fragmentation region. 
The resolved contributions decrease with p^ more fast, that is connected with 
the finite phase space value. 



6 Conclusion 

We have compared the /c^-factorization approach for heavy quark hadro- and 
photoproduction at collider energies with the existing experimental data. The 
agreement is reasonable when we use the cut-off (22) in Eqs. (13) and (14) and 
it does not practically depends on the value of QCD scale /x in these equations. 
We present also some predictions for the total cross sections of heavy quark 
production at LHC energies. 

Another example of very successful comparison of experimental data with 
the fc'T-factorization approach can be found in [10] where some different as- 
sumptions were used. 

It has been shown in Sect. 4 of [22] that the contribution of the domain 



with strong qt ordering (giT,?2T <^ rriT = Jit^q + Pt) coincides in the /ct- 
factorization approach with the LO PM prediction. Besides this a numerically 
large contribution appears at high energies in /c^^-factorization approach in 
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the region gi,2T > ""^t- This kinematically relates to the events where the 
transverse momentum of heavy quark Q is balanced not by the momentum of 
antiquark Q but by the momentum of the nearest gluon. 

An interesting consequence of such a kinematics is the fact that even at 
large transverse momentum piT » m the ratio of the single quark inclusive 
cross sections for charm and beauty production Rc/b = {dcfc/ (Pp) / {daf,/ (Pp) 
does not tend to unity. At Tevatron-LHC energies this ratio Rc/b — 1.5 and 
only selecting the events where both quarks has a large pt we obtain the ratio 
close to unit (the difference Rc/b — 1 does not exceed a few per cents when both 
PiT,P2T > 40 GeV). 

The configurations with qi^2T > ""^t are associated with the NLO (or even 
NNLO, if both Qit, q2T > ""^t) corrections in terms of the PM with fixed 
number of flavours, i.e. without the heavy quarks in the evolution. Indeed, 
as was mentioned in [1], up to 80% of the whole NLO cross section originates 
from the events where the heavy quark transverse momentum is balanced 
by the nearest gluon jet. Thus the large "NLO" contribution, especially at 
large pt, is explained by the fact that the virtuality of the t-channel (or u- 
channel) quark becomes small in the region q^ — pr, and the singularity of 
the quark propagator l/(p — q — iTig) in the "hard" QCD matrix element, 
M{qiT,q2T,PiT,P2T), reveals itself. 

The double logarithmic Sudakov-type form factor T in the definition of the 
unintegrated parton density Eq. (14) comprises an important part of the virtual 
loop NLO (with respect to the PM) corrections. Thus we demonstrate that 
the /cT-factorization approach collects already at the LO the major part of the 
contributions which play the role of the NLO (and even NNLO) corrections to 
the conventional PM. Therefore we hope that a higher order (in as) correction 
to the /cT-factorization could be rather small. 

Recall that our results are rather stable with respect of the factorization 
scale n variations and with respect to the form of angular cut-off parameter A 
in (22, (23). 

Another advantage of this approach is that a non-zero transverse momen- 
tum of QQ-sjstem {pxpair = Pit + P2T = qir + Q'2t) is naturally achieved in 
the /cj-factorization. We have calculated in [23, 24] the pxpair distribution and 
compared it with the single quark pt spectrum. At the low energies the typ- 
ical values of prpair are much lesser than the heavy quark px in accordance 
with collinear approximation. However for LHC energy both spectra become 
close to each other indicating that the transverse momentum of second heavy 
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quark is relatively small. The typical value of this momentum depends on 
the parton structure functions/densities. It increases with the initial energy 
and with the transverse momenta of the heavy quarks, pt- Thus one gets a 
possibility to describe a non-trivial azimuthal correlation without introducing 
a large " phenomenological" intrinsic transverse momentum of the partons. 

The contribution coming from qt > fi region could enhance the flux of 
colliding gluons and by this way partly explain the FNAL-Tevatron puzzle - 
data on the cross section of bb (or high-pj^ prompt photon) production are 2-3 
times larger than the conventional NLO PM QCD predictions [33, 43, 52] ^. 

At the moment we have no realistic unintegrated parton distributions which 
fit the data with accounting for the contribution from q^ > /i. We hope 
that future experiments will allow to distinct between different approaches. It 
is necessary to note is that the theoretical results concern the heavy quarks 
rather than the hadron production which can be investigated experimentally. 
The hadronization leads to several important effects, however their description 
needs additional phenomenological assumptions, see e.g. [53, 54, 55]. 

We are grateful to M.G.Ryskin for many useful discussions. This paper 
was supported in part by grants RCGSS-1124.2003.2 and NATO PDD (CP) 
PST.CLG 980287. 



^Recall, however, that the analysis [48] of the new CDF data gives a lower beauty pro- 
duction cross section. 
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Figure Captions 

Fig. 1. Low order QCD diagrams for lieavy quark production in pp (pp) 
collisions via gluon-gluon fusion. 

Fig. 2. Low order QCD diagrams for heavy quark production in direct 7p 
collisions via photon-gluon fusion. 

Fig. 3. The unintegrated gluon distributions integrated over q^ until 
Q^ = 100 GeV^ using Eq. (21) for A = fi/{fi + qr) (dashed curve) and 



A = nl J IJ? + q^ (dotted curve) with /x^ = Q"^ together with gluon struc- 
ture function GRV94 HO (solid curve) which was used for the calculations of 
these unintegrated distributions using Eq. (14) (a). The comparison of KMR, 
Eq. (14) (sohd curves), GLR, Eq. (11) (dotted curves) and KMS [39] (dash- 
dotted curves) unintegrated gluon distributions at x = 0.0005 (upper curves) 
and X = 0.05 (lower curves) (b). 

Fig. 4. The cross sections for beauty production for p^ higher than p™™ 
(a, c) and da/dE^ (b) in pp collisions at 1.8 TeV (a,b) with \yi\ < 1, and 0.63 
TeV \yi\ < 1.5, (c) and their description by the /c^-factorization approach with 
unintegrated gluon distribution fg{x, qr, jj) given by Eq. (14). The value of A 
in Eq. (13), (14) was taken equal to fJ,/{qT + jj) with the scale values /i^ = m^ 
(solid curves, jj? = 4m|, (dashed curves) and jj? = s (dotted curves). 

Fig. 5. The cross sections for charm production in pp collisions at 1.96 
TeV with |yi| < 1 and their description by the /c^^-factorization approach with 
unintegrated gluon distribution fg{x, qr, jj) given by Eq. (14). The value of A 
in Eq. (13), (14) was taken equal to fi/{qT + Ij) with the scale values /x^ = m\ 
(solid curves, /i^ = Am^ (dashed curves) and yU^ = s (dotted curves). 

Fig. 6. Total cross sections of charm (a) and beauty (b) photoproduction. 
The values of A in Eq. (13), (14) were taken equal to /i/(gr + A*) with scale 
values jj? = m^ (solid curves), jj? = 4m|. (dotted curves) and fi^ = s (dashed 
curves). The resolved photon contributions are shown by lower solid curves. 
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